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We prove that a beam splitter, one of the most common optical components, fulfills several classes 
of majorization relations, which govern the amount of quantum entanglement that it can generate. 
First, we show that the state resulting from k photons impinging on a beam splitter majorizes 
the corresponding state with any larger photon number k' > k, implying that the entanglement 
monotonically grows with k. Then, we examine parametric infinitesimal majorization relations as 
a function of the beam-splitter transmittance, and find that there exists a parameter region where 
majorization is again fulfilled, implying a monotonic increase of entanglement by moving towards a 
balanced beam splitter. We also identify regions with a majorization default, where the output states 
become incomparable. In this latter situation, we find examples where catalysis may nevertheless 
be used in order to recover majorization. The catalyst states can be as simple as a path-entangled 
single-photon state or a two-mode vacuum squeezed state. 
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I. INTRODUCTION 

In quantum optics, one of the most common trans- 
formations consists in the linear coupling between two 
modes of the electromagnetic field, as effected, for in- 
stance, by a beam splitter in bulk optics or an optical 
coupler in fiber optics [T] . Mathematically, it corresponds 
to a rotation in phase space, namely 

f a\ ( a'\ fcosO — sm0\ ( a\ 

[b)^{b') = { S me cose ){b)> (1) 

where a and b are the annihilation operators of the two 
modes that are coupled, while the angle G [0, 7r/2] 
is a coupling parameter related to the transmittance 
t = cos 2 9 of the beam splitter. The beam splitter is 
called balanced when t = 1/2 or = 7r/4. The trans- 
formation ([I]) belongs to the set of Gaussian unitaries 
as it corresponds, in the Heisenberg picture, to a linear 
canonical transformation in the annihilation (creation) 
operators and b^\ or equivalently to a quadratic 
Hamiltonian, namely H = i((vb — aw), see [2] for a re- 
view on Gaussian transformations. Moreover, it is a pas- 
sive linear operation since it conserves the total photon 
number, so it can be realized with a simple glass plate 
with a thin coating and no additional (pump) energy. 

The beam splitter is very conveniently modeled within 
the so-called symplectic formalism by focusing on the ac- 
tion of the rotation on the first- and second-order mo- 
ments of the quadrature operators. This enables treat- 
ing complex optical circuits made of beam splitters and 
other optical devices in a very concise way, which is suffi- 
cient for many purposes, e.g., when the goal is to predict 
distributions in phase space such as Wigner functions. 
However, if we want to make predictions about entropies 
or entanglement, we need to move back to state space 
and work with density operators. Such calculations are 
often nontrivial despite the simplicity of the transforma- 
tion in phase-space representation. For example, consider 



a single-photon state |1) impinging on a balanced beam 
splitter, the other input mode being in the vacuum state 
|0). The two-mode output state is obtained simply by 
inverting transformation ([I]) and writing the input-mode 
annihilation operators a and b as functions of the output- 
mode annihilation operators a' and V ' . The input state 
being a^|0), we can write the output state as 

2- 1 /2 ( a't + 6't)| ) = (|l)|0) + |0)|l))/^, (2) 

which is well known to be a path-entangled state of 
one photon, characterized by an entanglement entropy of 
1 bit as measured by the reduced von Neumann entropy 
of any output mode. However, whenever we consider 
higher photon-number states at the input and arbitrary 
transmittances, it becomes much harder to find closed 
formulas for the entanglement entropy. 

As another example illustrating the difficulty of treat- 
ing a beam splitter in state space, let us consider an arbi- 
trary input state in mode a that is coupled to a thermal 
field in mode b. The transformation a — > a' can be viewed 
as a thermal bosonic channel, which is a special case of 
a Gaussian phase- insensitive bosonic channel [3J. In or- 
der to derive the channel capacity, a crucial element is to 
determine the input state in mode a that results in the 
minimum-entropy output state in mode a! . Although it 
is very tempting to assume that this extremal input state 
is the vacuum state |0), this has not been proven to date 
[3]. It is linked to the Holevo- Werner conjecture, which 
states that Gaussian mixtures of Gaussian states achieve 
the capacity of such Gaussian channels [5]. 

In this paper, we exploit majorization theory in or- 
der to study the entanglement generated by an optical 
beam splitter. Majorization provides a partial order re- 
lation between bipartite pure quantum states and gives 
a necessary and sufficient condition for the existence of a 
deterministic LOCC (local operation and classical com- 
munication) transformation from one state to another 
O [7]- Here, we will show that a beam splitter obeys 
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two classes of majorization relations, which bear some 
similarity to those characterizing another optical compo- 
nent, namely a two-mode squeezer |S]. Specifically, we 
will prove that when feeding the input mode of a beam 
splitter with k photons while the other input mode is in 
the vacuum state, the resulting two-mode output state 
majorizes the state corresponding to any larger photon 
number k' > k. This implies that any bipartite entangle- 
ment measure on the output modes increases with k in a 
monotonic fashion. 

Then, we examine majorization relations with respect 
to the coupling parameter 9 (or, equivalently, the trans- 
mittance r) of the beam splitter. For a fixed photon- 
number input state \k) in one port and vacuum in the 
other port, we probe the existence of majorization re- 
lations between the output states corresponding to dif- 
ferent #'s, which we call parametric majorization. We 
find that in a region of finite width, the output state 
with parameter 9 majorizes the output state with 9' > 9, 
which implies that entanglement can only increase with 
the coupling between the two modes in this region. In- 
terestingly, we also disproof parametric majorization in 
other regions of the parameter 9, which means that the 
corresponding output states are then incomparable. In 
some cases, however, we can prove that these incompara- 
ble output states can be catalyzed [9], that is, if we sup- 
plement both states with an appropriate catalyst state, 
then the two product states become comparable (one is 
majorized by the other). Remarkably, the catalyst state 
can be as simple as a path-entangled single-photon state 
or a two-mode vacuum squeezed state. 

The rest of this paper is organized as follows. After 
summarizing the basics of majorization theory in Section 
II, we exploit it in Section III and prove that the input 
state |fc)|0) yields an output state |q/( fc )^ that majorizes 
the output state )) with k' > k. Then, in Section 
IV, we proceed to investigate parametric majorization, 
namely the relation between the output states \*f>( k \6)) 
and \ty( k '(9 + e)) corresponding to different transmit- 
tances. We focus, in particular, on infinitesimal ma- 
jorization relations, namely the case where e > is in- 
finitesimal. We observe a default of majorization in some 
specific regions of parameter 9, indicating that the above 
states are then incomparable. We also show how ma- 
jorization may be recovered by exploiting the concept 
of catalysis, namely by supplementing the state with an 
appropriate catalyst state. Finally, we close with the 
conclusions in Section V. 



II. THEORY OF MAJORIZATION 

The theory of majorization gives a means to compare 
two probability distributions and to conclude which of 
the two is more "disordered" or more "random" |10H13j . 
Consider two <i-dimensional real vectors p and q. We say 



that p is majorized by q, symbolized by p -< q, iff 

Erf < XX (3) 



i=i »=i 



for k = 1 , . . . , d — 1 and 



Erf = E< 



(4) 



where the down-pointing arrow on p and q indicates that 
the components are sorted in non-increasing order. Equa- 
tion Q is automatically satisfied if p and q are vectors 
normalized to unity, e.g., if they are probability distri- 
butions. Majorization only provides a partial order, in 
the sense that if p is not majorized by q (symbolized by 
p^q) then this does not imply that p >- q. When both 
p yi, q and q ^ p hold, we say that the two vectors are 
incomparable. 

The definition of majorization through Eqs. Q and 
Q is very handy for calculation purposes, but it does 
not makes it clear in what sense p is more disordered 
than q. A more intuitive definition is to say that p is 
majorized by q iff there exists a set of d— dimensional 
permutation matrices II„ and a probability distribution 
{t n } such that 



p = ^t„n„ 



(5) 



Roughly speaking, the above definition says that p is ma- 
jorized by q iff we can obtain p by randomly permuting 
the components of vector q and afterwards taking the 
average over all permutations. 

The equivalence between these two different definitions 
of majorization is implied by Theorems 1 and 2 below. 
Indeed, the notion of majorization is closely related to 
the notion of doubly stochastic matrices. A real d x d 
matrix D — [Dij] is doubly stochastic if all its entries are 
non-negative, and each row and each column sums up 
to 1. The following theorem gives the relation between 
majorization and doubly stochastic matrices. 

Theorem 1: p -< q iff there exists a doubly stochastic 
matrix D such that p = D • q. 

The set of doubly stochastic matrices of a given 
dimension is convex. All extremes points of this convex 
set are the permutation matrices n„, so any doubly 
stochastic matrix can be expressed as a convex combi- 
nation of permutation matrices. This is expressed in the 
following theorem. 

Theorem 2: The d x d doubly stochastic matrices D 
form a convex set whose extreme points are all the d x d 
permutation matrices n„. 

The convex set of d x d doubly stochastic matrices is 
called Bh'khoff's polytope. It admits d\ vertices (i.e., the 
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number of d x d permutation matrices) and its dimen- 
sion is (d — l) 2 . Note that if we want to express a point 
(a doubly stochastic matrix) belonging to this (d — l) 2 - 
dimensional polytope as a convex combination of the ex- 
tremal points, Caratheodory's theorem implies that we 
would need (d — l) 2 + 1 extremal points at most. 

One naturally expects that majorization theory should 
be connected with various measures of "disorder" , such 
as entropies. Indeed, since p -< q means that p is more 
disordered that q, any measure of disorder S : R d — s- K 
should satisfy 

S(p) > S(q) (6) 

for all p and q such that p -< q. A function S obeying this 
property is called Schur-concave. Consider, for example, 
the Shannon entropy 

d 

£i(p) = -X> ln ^ ( 7 ) 

or the Renyi entropy 

* Q(P)= l^ ln (5> Q ) (8) 

i=l 

of order a > 0, a ^ 1. (In the limit a — > 1, the Renyi 
entropy converges to the Shannon entropy.) These func- 
tions can be seen to be Schur-concave as a consequence 
of the following theorem [IT] . 

Theorem 3: p -< q iff J2i=i h (Pi) > J2i=i Kit) for au 
concave functions h. 

The usefulness of majorization in quantum informa- 
tion theory appears first if we wish to compare two 
density matrices instead of probability distributions. 
Consider the density matrices p and a of a d— level quan- 
tum system, and their respective vectors of eigenvalues 
A(p) and A(cr). We have the following theorem. 

Theorem 4: A(p) -< A(cr) iff state p can be obtained 
from state a by applying a mixture of unitarics. 

The proof goes simply by noting that there is a unitary 
transformation that aligns the eigenbasis of a with that 
of p, and that each permutation of the eigenstates can by 
realized by a unitary transformation. 

The connection with quantum information theory be- 
comes even more evident in the context of comparing 
pure bipartite entangled states. Indeed, majorization 
theory gives the means to determine whether one pure bi- 
partite state is convertible into another one using LOCC 
(local operation and classical communication). Consider 
two d— level systems A and B. Any bipartite pure states 
on these systems can be written in the Schmidt form 

l*)=E^>^' ( 9 ) 

i=l 



where and are suitable orthonormal bases of 

the systems A and B, respectively. The reduced density 
matrix of system A is = trs|^')(^| = Yli=i ^N)(*U, 
and similarly for B (the two reduced density matrices 
have the same eigenvalues A^). We have the following 
theorem [SJ [7]- 

Theorem 5: State \^) can be converted determinis- 
tically into state |$) by means of LOCC iff A* -< A$, 
where A^, is the vector of eigenvalues of p\ = trs|^)(^ , | 
and similarly for A$. 

For conciseness, we will write simply \E' -< <E> instead of 
A^ -< A$ in what follows. A consequence of this theorem 
is that \fr -< $ iff > p(&) for all measures of entan- 
glement p. A measure of entanglement, or entanglement 
monotone, is a non-negative function of the state which 
does not increase on average under LOCC and vanishes 
on separable states [14] . A common example is the en- 
tropy of entanglement, E(ty) = Si (A*). Since, accord- 
ing to Theorem 5, converting into $ is possible when 
* -< $, this means that p(^) > p(&) must hold for 
all entanglement monotones. Conversely, the maximum 
probability of success of converting \I> into $ by means 
of an LOCC protocol satisfies P(* -> $) < min p 
where the minimum is taken over all entanglement mono- 
tones [T3]. Since a strategy exists where P(\& — > <I>) = 1, 
this implies that p(^>) > /«($), V/x. 

According to Theorem 5, if A* and A<j, are incompa- 
rable, then there does not exist a strategy to convert one 
state into the other by LOCC with probability 1. Re- 
markably, it has been shown that one may nevertheless 
be able to accomplish such a transformation dctcrminis- 
tically with the use of an auxiliary entangled state, an 
effect called catalysis [5]. If two states |W) and |$) have 
incomparable A vectors, then, under certain conditions, 
there exists an entangled state |C) that the two parties 
can share, called a catalyst state, such that 

l*> ® \c) ^ c 1$) ® \C) (10) 

is possible. The term "catalysis" is justified because the 
entangled state |C) is recovered after the LOCC trans- 
formation. Note that if converting iff into <& by catalysis 
is possible, then all additive measures of entanglement 
must satisfy /x(\&) > p(<&)- In particular, we must have 
£a(A*) > S a (\J) for all a > 0. 

III. MAJORIZATION WITH RESPECT TO 
PHOTON NUMBER 

We now prove a first class of majorization relations 
characteristic of a beam splitter by using the definition 
of majorization involving doubly stochastic matrices. Let 
\iff( k \6)} be the output state of a beam splitter if the 
input state is |fe,0), as shown in Fig. 1. Denoting by 
U(8) the unitary transformation resulting from the beam 
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|*M(0)) 



lfc+1) 




This simply expresses the fact that the probability of 
transmitting n photons out of k + 1 incident photons is 
the sum of two mutually exclusive possibilities: either 
n — 1 photons (out of k) are transmitted and the (k + 
l)-th photon goes through, or n photons (out of k) are 
transmitted and the (k + l)-th photon is reflected. We 
can expand Eq. ( 15 ) as 



P 



>(*+!)/ 



(9) = + 



P^ +L> (9) = cos 2 9P[ k) {9)+s\n 2 9P [ 2 k) {9), 



P^ L> [6) 
,(fe+i 



FIG. 1: Majorization relations with respect to the input pho- 
ton number. In the upper scheme, the input state is |fc,0), 
while in the lower scheme it is \k + 1,0). The output state 
|* (fc) (6»)) majorizes |* (fc+1) (6»)}, so that the generated entan- 
glement n(ty( k ' (9)) monotonically grows with k for all 6. 



splitter, we have 



- U{9) \k,0) = Y, \/Pn\8) \n, k-n), (11) 

n=0 

where 



Pi k \9) = 



cos 



2n a sin 2(fe-") Q_ 



(12) 



The reduced density matrix of the first output mode is 

k 



p«(0) = £p«(0)|rz)H 



(13) 



where P^ ] (0) can be interpreted as the probability that 
n photons are transmitted out of the k incident photons 
if the transmittance of the beam splitter is r — cos 2 9. 
We wish to prove a majorization relation between 

Prt\9) and Pn (0), that is, we want to prove that 
there exists a doubly stochastic matrix D such that 

p(*+i)(0) = D (fc+1) • P (fe) (6>), (14) 

where P (fe) (6») is a vector having the eigenvalues P^ (9) 
of // fe ) (9) as components. Note that in Eq. (14), we wrote 



p( fe+1 ) on the left-hand side and on the right-hand 
side, which means that p( fc+1 ) is more disordered than 
P^ k \ Using Pascal identity for the binomial coefficients, 
we obtain the recurrence equation 



Pi k+1) (9) = { * ^ 1 ) cos 2 " 9 sin 2 ( fc+1 -«) , 




cos 2 " sin W""), 



P^(B) cos 2 8 



(15) 



p&YHe) 



cos 2 9P { k k) (9)+0. 



(16) 



This can be put in the form of Eq. ( 14 ) by defining 



>(fe+i) 



/P o (fe+1) (0)\ 
Pi fe+1 )(0) 



Pf +1) (0) 



P (fe) (0) 



/P o (fc) (0)\ 
P 2 (fe) W 



pf (0) 

V o ) 



,(17) 



where a zero entry has been inserted in the vector pw in 
order to make and p( fc+1 ) of equal dimension. The 
doubly stochastic matrix is 



D (fe+i) 



/ sin 

cos 










2a sin 2 9 
cos 2 sin 2 



V o 











cos 2 8\ 





(18) 



sm' 



where the right-most entry of the fist row has been chosen 
so as to fulfill the doubly-stochastic conditions (it plays 
no role since the last entry of P ( - fe - ) is zero). 

Thus, we have proven the majorization relation 

<Er(*+l)(0) ^ 4,(fc)(0) 5 



V0, 



(19) 



which implies that when increasing the number of the 
incident photons, the 2-mode output state can only be 
more entangled and the corresponding 1-mode reduced 
states are getting more disordered. This also implies 
that the two parties (Alice and Bob) can convert state 
|\j/( fc +!)(0)) i n to state |*W(0)) by using a deterministic 
LOCC transformation. To achieve such a transformation, 
Alice can, for example, perform a two-outcome POVM 
measurement with the following Kraus operators, 



^ = E 



n=0 



7 (fe) jk + l-n 



k+1 



n)(n\ 



and 



(20) 



(21) 
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satisfying J-\j-\ + T\T-z = 1 for all k. Then, she must 
communicate her outcome to Bob, who has to apply 
proper local unitaries. If outcome 1 occurs, then Bob 
should apply the unitary 

k 

U[ k) =5^|n)<n + l| + |fc + l)(0| (22) 
n=0 

corresponding to a cyclic shift in the space spanned by 
{|0), • • • \ k + 1)}. The second term on the right-hand side 
of Eq. p2| ) ensures unitarity (it plays no role since Bob's 
reduced state is supported by {|1), • • • \ k + 1)} when out- 
come 1 occurs). If outcome 2 occurs, then Bob does noth- 

(k) 

ing, that is, he applies the unitary = U. It is easy 
to check that the transformation |*( fc+1 )(0)) ->■ |*( fc )(6»)) 
takes place for both outcomes, so the LOCC transforma- 
tion is indeed deterministic. 




FIG. 2: Majorization relations with respect to the coupling 
parameter 9 (or transmittance r = cos 2 6). In both schemes 
the input state is \k, 0) but the angles differ by e. In some 
parameter region, the corresponding output states |\I/' fc '(0)) 
and \^/^ k \9 + e)) are proven to satisfy a majorization relation. 



IV. PARAMETRIC MAJORIZATION WITH 
RESPECT TO TRANSMITTANCE 

A. Infinitesimal majorization 

We now examine the scenario that is summarized in 
Fig. 2. The input state is fixed to |fe,0), but we change 
the angle 9 parameterizing the transmittance of the beam 
splitter by an infinitesimal amount e. Note that we take 
> 0, e > 0, and 8 + e < j. (For angles greater than ~, 
the transmittance and reflectance just interchange their 
roles.) An equivalent way to see this scenario is depicted 
in Fig. 3. Our goal is thus to probe whether the inter- 
mediate state \^ k \0)) majorizes or not the final state 
\VW(6 + e)). To this end, we find it easier to use the first 
definition of majorization, involving the accumulations of 
the ordered vectors of eigenvalues of the reduced density 
matrix. We will refer to these vectors as OSC (ordered 
Schmidt coefficients). 

Let P^(8) be an OSC vector, whose components are 
the elements of the binomial distribution ( 12 ). From now 



on, we drop the index k as it is fixed. This OSC vector 
will not keep the same ordering as the parameter 6 varies, 
so we will adopt the notation P^ r (8), where r — 1,2,... 
labels the regions of parameter 8 in which the ordering 
of the OSC vector remains the same. More precisely, 
we have a change of ordering every time two eigenvalues 
P n (8) and P m (8) are equal, which occurs at 



= arctan 



(k — n)!n! 
[k — m)!m! 



(23) 



Indeed, it can be shown that if P n {9) = P m {9) for m ^ 
n, then dP ^g S ^ ^ > i- e -7 t ne two eigenvalues cross. 

The cross-over points between regions are the different 
solutions 8\ < 82 < ... of the above equations. We define 
the region r — 1 as 8 e [0, the region r = 2 as 
6 [81,82), etc. 



I*>- 



10} 





|tf«(0)) 



|* (fe) (0 + e)> 



FIG. 3: Same situation as in Fig. 2, viewed as a sequence of 
two beam splitters with angles 9 and e. In some parameter 
region, \¥ k) {0)) majorizes \¥ k) (9 + e)) for all k. 



Our goal now is to check whether the infinitesimal ma- 
jorization relation 



P ir (8 + e) -< P ir {8) 



Irr 



(24) 



holds or not within region r, taking the limit of an in- 
finitesimal angle e. Using the definition of Eq. ([3]), we 
have to prove 



E P n r ( d + £ ) < E P n r (0), 3 = 0, ■ • ■ k - 1 

(25) 



n=0 n=0 

^E^^°> j = o,...*-i 



n=0 



dd 



By defining the vector of accumulation derivatives 

' dP^(8) 



- E 



d0 ' 



(26) 



the infinitesimal majorization relations can be written 
simply as 



4 r (0)<O, j = 0,...k-l. 



(27) 



We do not need to consider the last accumulation deriva- 
tive al r (9) = since Eq. Hh is necessarily verified (the 
OSC vectors are normalized). 



G 



The violation of at least one relation in Eqs. (27) is 



sufficient to disproof majorization in region r. A priori, if 
the above majorization relations do not hold, there may 
nevertheless be a majorization in the opposite direction 
if all relations are satisfied with > instead of <. However, 
the (k — l)-th accumulation is the same in all regions no 
matter what the ordering is, and its derivative 



a k -i(0) = -2ksm 2k ~ 1 6 cos 



respects Eq. (27) with a strict inequality (except in the 
trivial cases k — or 9 = 0). Hence, majorization is 
never possible in the opposite direction. 

It is easy to see that in the region r = 1, the compo- 
nents of the OSC vector are 



Pf(0) = 



sin 2n 6»cos 2(fc - 



(29) 



This region extends until the second largest eigenvalue 
becomes equal to the largest eigenvalue and the two 
switch places. Beyond this cross-over point, we enter the 
second region r = 2. Now, let us prove that the infinites- 
imal majorization relation P^ 1 (9 + e) -< P^ 1 (0) always 
holds in region r = 1. We have 



de 



2n 
tan£ 



2{k — n) tan( 



and 



Pk\e) = Ftf{e) 



tan 



2" , 



Pi 1 (6) (30) 



(31) 



from which we can express the vector of accumulation 
derivatives as 



a i (0) J ' 1 = P O U (0) [ 2n ~ 2 ( k ~ n ) tan2 9 

re=Q 

k 



tan 



2n-l 



(32) 



The summation in Eq. ( 32 ) can be expressed in a closed 
form as 



a 3 {9) 



1 1 - -Pf(0) 2(k - j) r ) tan 2 ^ 1 0, (33) 



which is non-positive for j = 0, . . . k — 1. Thus, infinitesi- 
mal majorization holds within region r = 1, which means 
that all states are comparable within this region, 



* (fe) (0 + e) -< # (fe) (0), Vfc>0, 



(34) 



and all measures of entanglement increase with 9. 

At the cross-over point between regions r = 1 and 
r = 2, the first two components of the OSC vector switch 
places and the derivative of the first accumulation be- 
comes dP ^g e ^ , which is positive at this point. Therefore, 
majorization is violated from the left boundary of region 



r = 2 until the point where this derivative ceases to be 
positive (and possibly beyond that point). In general, 
in every region that begins with a positive derivative of 
the first accumulation, which is equal to for some 

n, majorization is violated at least until the point where 
this derivative ceases to be positive. From Eqs. (30) and 



(|31|), we find that the derivative dP 2g 6 ^ remains positive 



(28) up to the value 9+ < arctan i 



In order to illustrate how parametric majorization be- 
haves, let us exhibit three examples, corresponding re- 
spectively to a single-photon, two-photon, and three- 
photon state impinging on the beam splitter. 



B. Examples 

Example 1: 

We first consider the case of a single photon (k — 1). The 
OSC vector in region r = 1 is 



(35) 



In order to find all possible regions we have to find all 
solutions of cos 2 9 = sin 2 9 in [0, j). The only solution is 
9 = ~, which means that there is a single region r = 1 
and, as we proved earlier, parametric majorization holds 
everywhere. 

Example 2: 

We now move to the case of two photons {k = 2). The 
OSC vector in region r = 1 is 



cos t) 
2 cos 2 9 sin 2 9 
sin 4 9 



(36) 



where this ordering holds for [0, arctan ^). There is 
a second region r — 2 corresponding to [arctan , f ) , 



where the OSC vector is 



'2 cos 2 6 sin 2 N 
P 4 - 2 (6») = I cos 4 6> 



(37) 



As proven in full generality, majorization holds in re- 
gion r = 1. However, in region r — 2, the accumulation 
derivatives are given by 



af(9) = 2 cos* sin* 
af(0) = -4 cos 9 sin 3 i 



(38) 



The accumulation derivative a (9) is positive in [0, J), 
so majorization does not hold in the entire region r = 2. 
This means that there ought to be measures of disorder 
that decrease instead of increase as a function of in re- 
gion r = 2. Indeed, we observe in Fig. 4 that although 
the Shannon entropy increases, all other Renyi entropies 
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of order a > 1 exhibit a decreasing behavior somewhere 
within the region r — 2. In particular, the Renyi en- 
tropy of order a — > oo, which is the min-entropy and 
is directly related to the leading probability of the OSC 
vector, starts decreasing immediately when we enter the 
second region at = arctan . 




FIG. 4: Entanglement Renyi entropies resulting from a 2- 
photon state impinging on a beam splitter as a function 
of the coupling parameter 8 (related to the transmittance 
t — cos 2 8). The vertical dashed line denotes the bound- 
ary between parameter regions r = 1 and r = 2. The von 
Neumann entropy (a — > 1) keeps increasing in region r = 2, 
while higher-order Renyi entropies have a different behavior 
and start decreasing somewhere in region r = 2. The min- 
entropy (a — > oo) exhibits a non-differentiable point right at 
the crossover point and decreases throughout the entire region 
r = 2, reflecting the default of majorization. 

Example 3: 

As a last example, we consider the case of three photons 
(k = 3). We have two cross-over angles, 



0i 
2 



arctan 



arctan 



V3 
1 

w 



(39) 



which define three regions of different orderings in [0, J) 
and three corresponding OSC vectors, 



I cos 6 \ 
3 cos 4 9 sin 2 9 
3 cos 2 9 sin 4 6> 
sin 6 



P i2 (0) 



P i3 (0) 



/3 cos 4 sin 2 0\ 
cos 6 
3 cos 2 sin 4 
sin 6 

/3 cos 4 sin 2 
3 cos 2 sin 4 



V 



cos" 
sin 6 



(40) 



In region r = 1, it is easy to confirm that majorization 
holds, as it should. In regions r — 2 and r = 3, the 



accumulation derivatives are given, respectively, by 
af{9) = 3cos 3 0(-l + 3cos20)sin0 



and 



af{9) 
af{9) 

af{9) 
af{9) 
af{9) 
The quantity 



= --sin 3 20 
= -6 cos sin 5 



(41) 



3cos 3 0(-l + 3cos20) sin t 
3 

- sin 40 
2 

-6 cos sin 5 0. 



(42) 



a 2 



0, arctan 



V2> 



(0) is positive in the interval 
which means that we are sure that there 

is no majorization in the interval [arctan -^=, arctan 
i.e., from the left boundary of region r — 2 up to where 
<2q 2 (0) remains positive. In region r = 3, a^{6) is always 
positive, while the other accumulation derivatives are 
negative within this region. Hence, for r — 3, the 
states are always incomparable. This last violation of 
majorization is, however, not visible with the Renyi 
entropies. In Fig. 5, we display the evolution of entropies 
across the three regions. 




FIG. 5: Entanglement Renyi entropies resulting from a 3- 
photon state impinging on a beam splitter as a function of 8. 
The two vertical dashed lines at 8\ and 82 separate the three 
regions r = 1,2,3, while the dotted line corresponds to the 
local minimum of the min-entropy. Majorization is violated in 
the region r = 2 from the left boundary of this region at 82 up 
to the local minimum of the min-entropy. The majorization 
violation throughout the entire region r = 3 is not manifested 
by the behavior of the Renyi entropies. 



C. Catalysis 

We have shown that one can always expect a default 
of majorization when the leading term in the probability 
vector changes, and this majorization default prevails at 
least up to the point where the first accumulation deriva- 
tive ceases to be positive, or, equivalently, until we reach 
the local minimum of the min-entropy within this region. 
Beyond the case where the first two components of the 
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Y Y 

*< 3 >(0.62)) 18 |C(0.7)) |*< 3 )(0.72)> g> |C(0.7)> 

FIG. 6: Schematic of the catalyzed conversion between the 
incomparable states resulting from a 3-photon state impinging 
on a beam splitter with angles 6 — 0.62 and 6 + uj = 0.72. The 
catalyst is the entangled state obtained from a beam splitter 
with angle 6 = 0.7 and one single-photon input state. 



probability vector switch places, it appears difficult to 
provide general rules for predicting the existence or ab- 
sence of majorization for an arbitrary k, and one has to 
treat the problem on a case-by-case basis. The situa- 
tion also becomes more complicated if we take a non- 
infinitesimal angle e such that the pair of angles 9 and 
9 + e belong to different regions. 

It is natural to ask whether the incomparable states 
that occur when we change the parameter 9 can never- 
theless be made comparable through catalysis. We will 
show that this is indeed possible in certain cases, and will 
provide an example for this. Note that not all incompa- 
rable states can be catalyzed: some necessary conditions 
have to be respected [3 fTSHTT] . To solve the problem of 
whether catalysis is possible or not in generality is dif- 
ficult due to the fact that one has to reorder the vector 
resulting from the tensor product of the state to be cat- 
alyzed and the catalyst state. 

Consider two angles 9 and 9 + lo that give OSC vectors 
P^ r (6*) and (9 + uj) which are incomparable. These 
angles may be within different ordering regions, which is 
the case in the following example where r = 2 and r' = 3. 
We take k = 3, 9 = 0.62, and ui = 0.10, which gives 



P i2 (0.62) 



and 



P i3 (0.72) = 



0.44439 
0.290641 
0.226491 
vO.0384782/ 



0.416698 
0.320544 
0.180565 
vO.0821927/ 



(43) 



(44) 



such that P ;3 (0.72) ^ P ;2 (0.62). The path-entangled 
single-photon state cos#|l)|0) + sin6*|0)|l) resulting from 
a single photon impinging on a beam splitter with angle 
9 = 0.7 is sufficient to serve as a catalyst for these two 
states. It corresponds to the binary probability vector 



C(0.7) 



0.584984\ 
0.415016 J ' 



(45) 



and one can easily verify from Eqs. (43), (44 1, and (45 1, 
that P i3 (0.72) <g) C(0.7) ~< P i2 (0.62) <g) C(0.7), implying 
that the catalyzed conversion is possible. This is sum- 
marized in Fig. 6. We could also use as a catalyst a 
two-mode squeezed vacuum state oc X^^Lo t an h n r \ n )\ n ) 
with squeezing parameter r = 1.38. Several other numer- 
ical examples can be found and some of them, like the 
ones provided above, are experimentally accessible. 

Note that, due to the additivity of the Renyi entropies, 
one should look for catalyzable incomparable states only 
in regions where all of the Reniy entropies increase (see 
also Refs. [H3III])- However, in the case at hand, we have 
got numerical evidence that the sole behavior of the min- 
entropy seems to give a necessary and sufficient condition 
for the existence of catalysis. The latter property will be 
examined in a forthcoming work. 



V. CONCLUSION 

We have found several classes of majorization relations 
characterizing a beam splitter, or more generally the lin- 
ear coupling between a pair of bosonic modes. More for- 
mally, we have proven that the passive Bogoliubov trans- 
formation of Eq. ([!]) fulfills some majorization relations, 
which enable comparing the output states correspond- 
ing to various input photon numbers k as well as various 
coupling parameters 9 (or transmittances r = cos 2 9). In- 
terestingly, this behavior is reminiscent of the majoriza- 
tion relations that have recently been shown to prevail 
with an active Bogoliubov transformation (a two-mode 
squeezer or parametric amplifier) 8 . Note that in con- 
trast to Ref. jS], the present analysis has no implications 
on the Holevo- Werner conjecture. Indeed, we consider 
situations where one of the input modes is in the vac- 
uum state, corresponding, in the language of quantum 
channels, to the case of a pure lossy channel. For that 
specific channel, the conjecture is well known to be true 
|20) since the vacuum input state results in another vac- 
uum state at the output of the channel (the vacuum is 
the extremal input state since the corresponding output 
entropy is minimum — it is zero). 

We have shown that for any value of the transmittance 
parameter 9, the output states resulting from injecting 
Fock states |fc) in one port of the beam splitter and vac- 
uum |0) in the other port obey a chain of majorization 
relations *( fe+1 )(6») -<! tf<W(0), for all k > 0. As a con- 
sequence, the output states can only be more entangled 
when increasing the number of incident photons, and we 
have found an explicit deterministic LOCC transforma- 
tion that maps *( fc+1 )(0) onto #( fc )(0). 

In contrast, we have found that the situation is more 
complicated when varying the parameter 9 and keeping 
k constant. In that case, we have shown that there exists 
a first region in the space of parameter 9 where a para- 
metric infinitesimal majoriation relation holds, by taking 
the limit of an infinitesimal angle e > for any k > 0, 
namely <]/( fe )(# + e) -< ty( k \9). This implies a monotonic 
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increase of the entanglement of the output states when 
decreasing the transmittance and moving towards a bal- 
anced beam splitter. However, beyond some value of the 
parameter 9, we have shown the existence of a default 
of majorization, which occurs because the ordering of 
the OSC vectors changes in such a way that the leading 
probability is replaced by another one. This majoriza- 
tion default holds from the left boundary of this new 
ordering region at least up to the local minimum of the 
min-entropy. Moreover, by examining specific examples, 
we have shown that one may find more violations of ma- 
jorization for non-infinitesimal angles e within the same 
ordering region or between different ordering regions. 

Finally, we have provided an example of two incompa- 
rable states, resulting from different values of 9, whose 
conversion can nevertheless be catalyzed with the help 
of an experimentally accessible state, such as a single- 
photon path-entangled state or a two-mode squeezed vac- 
uum state. Catalysis schemes like the one in Fig. 6 may 
potentially be used for authentication protocols based on 
entanglement-assisted LOCC [2T], [22] . Further investiga- 



tions should also include a more general solution to the 
catalysis process in the parameter-varying case, the anal- 
ysis of majorization relations in more complicated opti- 
cal circuits in the spirit of |23j . or more ambitiously the 
study of phase transitions and critical phenomena in a 
field-theoretical approach |24j under the prism of para- 
metric majorization, where the parameter could be the 
temperature of a thermal field 25 . 
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